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Bernard Kaufman 
ABSTRACT 

The perturbative effects of the four largest moons of Jupiter on a space- 
craft that passes close to Jupiter are investigated using a total of twenty dif- 
ferent encounter trajectories with the planet. An accurate model for the moons 
is utilized in the equations of motion which are numerically integrated over the 
time that the probe spends within Jupiter’s sphere of influence. It is shown 
that the perturbations caused by the moons do have effects on the post encounter 
heliocentric trajectories which are large enough to make necessary the inclu- 
sion of the moons in any detailed study to be undertaken. On the other hand, 
actual mission objectives such as the time to 10 astronomical units seem to be 
little affected by the presence of the moons. 

A valuable by-product of this study has been the programming of an accurate 
ephemeris for the Galilean moons which has not previously been available in this 
country as far as the author has been able to determine. 
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INTRODUCTION 


The ultimate goals of deep-space investigations are the understanding of 
the origin, laws and evolution of our solar system and of the galactic medium 
itself. By means of a probe to Jupiter and beyond into deep space we may ex- 
tend our investigations many fold. Such a probe is currently being planned by 
NASA and will be entering a region never before investigated which will include 
the asteroid belt and the trapped radiation region at Jupiter. 

A probe that is to be sent into deep space requires a tremendous amount of 
energy which with today's boosters can only be accomplished at the expense of 
the final payload. An imaginative way around this problem is to use the gravi- 
tational attraction of a planet. When a probe is in a free fall trajectory about 
the sun in such a manner that it passes within the vicinity of a planet, the tra- 
jectory is considerably changed by this encounter. The changes must, of 
course, depend on the encounter geometry of the probe and the planet. 

This report considers a probe that will utilize a flyby of the planet Jupiter 
to obtain changes in the post encounter heliocentric trajectory. A by-product of 
an encounter with Jupiter is the opportunity to study the planet itself. 

Many authors have treated a Jupiter flyby in the past (for example see 
references 8 thru 13) based on two body motion, however, Jupiter has twelve 
natural satellites, four of which are rather large and their masses are rela- 
tively well known. To the author's knowledge no one has done an integrated 
trajectory study which includes perturbations by the moons of Jupiter. In this 
report only the four Galilean moons are included since these are the largest and 
their masses are well known. The problem of obtaining accurate positions of these 
moons has resulted in the programming of Sampson's method to generate an ephem- 
eris of the moons. Such a program has not previously been available in this coun- 
try as far as can be determined. 
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THE TRANSFER TRAJECTORY 


The first problem that must be considered is that of the transfer trajectory 
from the earth to Jupiter. This will involve the determination of the proper 
launch conditions, parking orbit duration and injection energy to satisfy total 
flight time and nominal targeting at Jupiter. 

Many authors have treated this problem using various methods to obtain a 
solution. References 1 through 7 in the bibliography give some of these methods 
and also indicate the complexity of the problem. 

Since in the present study we are mainly interested in the probe once it is 
in the vicinity of Jupiter, the transfer orbit is only of minor interest. For this 
reason many of the constraints on the trajectory were dropped in order to 
simplify the solution. An approximate transfer trajectory is obtained using 
only the flight time and the positions of the launch and target bodies as con- 
straints. No attempt is made to refine the approximate trajectory in any way. 
The solution of the transfer trajectory under these restrictions is then solved 
using the well known solution to Lambert's problem. At the target planet the 
trajectory is altered by using the patched conic assumption which states that at 
any given time the probe is affected by one and only one body. At any time the 
gravitational field of a planet is considered to have a force field called the 
sphere of influence whose radius s is given by (reference 1) 



where m and M are the masses of the planet and sun respectively and R p is the 
planet's distance from the sun. 
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During the transfer phase of the probe, the Earth and Jupiter are consid- 
ered to be massless points until the probe comes within the sphere of influence 
of Jupiter. The method used to define this conic transfer trajectory is es- 
sentially that found in reference 7 and is described below. 

Solution of the problem of a transfer trajectory necessitates knowing the 
heliocentric position vectors of the Earth ( r* e ) at launch date and of Jupiter 
( 7 . ) at arrival date. These vectors can be obtained from an ephemeris or by 
using two body equations and a given position of the planets at some epoch. 
Figure 1 shows the essential transfer geometry. 



Figure 1 -Launch Geometry 

The heliocentric transfer angle 'P in figure 1 can be calculated from 7 e and 
7 . and is the angle between the position vectors measured in the direction of 
planetary motion. Two other quantities in figure 1 to be defined are the angles 
6 and y . Theta ( 6) is the true anomaly of the Earth at launch and y is the 
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heliocentric flight path angle which is the angle between a perpendicular to 7 e 
and the velocity vector v of the probe at launch. The angle 7 is the parameter 
that we wish to calculate. 

In order that the launch body and target body lie on a single conic orbit it 
is necessary that the following be true 

r = P (1) 

e 1 + e cos 8 

r . = P (2) 

J 1 + e cos (<9 + *P) 

where /0 = a(l-e 2 )is the semi latus rectum. 

From (2) 


p = r . [l + e cos 0 cos 'P- e sin 0 sin V] 


and frpm ( 1 ) 


e cos 8 - 



From figure (1) 

7 • v 

Tan 7 = cot (90 - 7) = — — ^ 

I r e x v | 


r • v 

e 

Va(l - e 2 ) 


_ e sin E _ r e e sin 6 _ 

P 

Therefore 


e 


s i n 9 - JL- Tan y 
r 

e 


substituting (4) and (5) into (3) we obtain 

r e r. (1 - cos V) 

p — ^ 

r — r. cos^+r. sin f Tan 7 
e j j / 


(3) 


(4) 


( 5 ) 


( 6 ) 
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Equation 6 characterizes the conic section connecting the two radius vectors 
and it is obvious that p is a function only of the unknown angle y . 

The calculation of y is an interative procedure and only one value of y will 
yield the desired flight time between the two planets. Reference to figure 2 
shows that two transfer trajectories of very short time exist. One is the 
straight line solution and the other bends very close around the sun (hot line) . 

A solution lying inside the triangle formed by these two is forbidden since no 
conic section with the sun at a focus exists in this region. Thus a lower limit 
(7 m i n ) can be calculated for y by using the straight line solution and it is easily 
seen that 


y . - Tan 

7 min 


-1 


*r . cos 'P - r 
J e 


n ^ 


(7) 


j j 

The second limit for y is that for a posigrade parabolic orbit for an infinite 
transfer time. Reference to figure 3 shows the typical geometry for this orbit. 



Figure 2-First Limit for Lounch Trajectories 



V 



Figure 3-Second Limit for Launch Trajectories 


Recall equation (5) 


e sin 6 - Tan 7 
r 


but for this orbit e = 1 and 


P =■ 


2 r 


1 + Tan 2 — 


Therefore 


. 6 6 

sin — cos — = 
2 2 


Tan 7 


2 e . 2 e 

co s ^ — + sin z — 

2 2 


cos 2 — 


_ 2 @ t 

- cos z — Tan 7 


or 


Tan — = Tan 7 

2 


Therefore for a parabola 


0 = 2 y 


From equations (1) and (2) we have 

r e (1 + cos 6) = r . [1 + cos (6 + 'F)] 


or 


r e - = ( r j cos'P-r e )cos<9-(r. sinV)sini9 

where & is to be eliminated. 

The square root of the sum of the squares of the coefficients of 0 is 

(r? + r 2 - 2r e rj cos *F) 1/2 = chord 


thus 


r - r . 

e — = cos a cos 6 - sin a sin 0 


chord 


where 


Thus 


r . cos W - r 

j e 

COS a = 

chord 




> 


sin a = 


r . sin V 
j 


chord J 


r - r . 

— 1 = cos (a + 6) - cos /3 

chord 


and cl can be evaluated from eqn (9) 
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Therefore 

and from (8) 
Thus 


a + 8 - (3 

Q - 2 y. 

I 3 - a 

^LIM = 2 


( 10 ) 


Equations (7) and (10) now allow us to compute the transfer trajectory in 
the following manner:* 
the initial guess for y is 


(y . if ¥ < 180° 

7 min 

I-tt/2 if ¥ > 180° 


and a step size (ss) for y. is chosen such that 

= Vl + SS 

and 


From eqn (6) 


P= 


y . < y < y . . 

7 mm 7 n 7 lim 


r e r. ( 1 -cos ¥) 


r . s in ¥ Tan y - r . cos ¥ + r 

i ' n j 


The velocity is given by 


2 2 Ms M s 
= 

r e 3 

PP S 


r 2 cos 2 y 

e 7 n 


where m s is the gravitational constant of the sun. 


*These two limits on y also rule out retrograde transfer orbits 
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Define 

C3= = .'!» 

r . 

which is twice the total energy per unit mass. 

It can easily be shown that 

e sin ^ = vsin y n 

which yields 6 . 

Then 



Appendix A shows how the time ( t e ) from periapsis to a true anomaly of 
6 can be calculated. In a similar manner the time ( t .) from periapsis to 
Jupiter can be calculated. Then the flight time is 

At = tj - t e . 

If At is within a predetermined tolerance to the desired flight time then 
nothing more need be done. If not, the method of false position allows an 
iterative solution to y for the desired flight time. The quantities finally ob- 
tained are v, a, e, 6 and y . 

Using the above calculated y we may calculate the velocity vector v by 
decomposing onto two unit vectors r® and (r e xr. )° x ? e ° as in figure 4. 

We obtain 

v=f [cosy (r e x fj)° x r e + sinyTJ 



( 11 ) 
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Figure 4-Launch Velocity Diagram 


Equation 11 represents the velocity vector of the probe when it is launched 
from the Earth towards Jupiter. This velocity vector along with W, the helio- 
centric transfer angle, and the position vector of the probe is all that is neces- 
sary to find the position and velocity of the probe at arrival at Jupiter. It is 
easily seen that for this approximate transfer trajectory, the position vector of 
the probe at launch is that of the Earth which is obtainable from the ephemeris. 
Appendix B describes the method used to update the position and velocity of the 
probe to the arrival date at Jupiter. 



UNDISTURBED JUPITER FLYBY AND POST ENCOUNTER TRAJECTORY 

The Sphere of Influence 


As discussed in the last section the heliocentric velocity vector of the 
transfer orbit can be obtained at the arrival date at Jupiter. Since in this ap- 
proximate scheme Jupiter is considered to be a massless point this arrival 
date is actually an impact time. The problem now is to determine the position 
and velocity of the probe at arrival at the sphere of influence of Jupiter. The 
heliocentric position vector of the probe at impact is obviously that of Jupiter. 

Using the equation for the radius of the sphere of influence ( S ) discussed 
earlier 


S = 



R 


p 


the method of false position and two body equations allows us to back up the r 
and v from impact to a distance S and thus obtain r and v at entrance to the 
sphere of influence. 

We let r and v be the heliocentric impact position and velocity. From r 
and v we can easily compute the semi-major axis a , the eccentricity e , the 
mean motion n and the true anomaly 0. at impact. From 6*. we find the 
eccentric anomaly E from 


„ yl-e^sinif?. 
Sin E = I i. 

1 + e cos 6 i 


Cos E = 


cos 9 i + e 
1 + e cos 0 . 


> 




( 12 ) 
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Kepler’s equation then allows us to compute the time t i from periapsis to < 9 . 


M = E-e sin E = nt. 

The next step is to start the false position interation by choosing a time 

t = t A -At n 

where At can be any arbitrary time. 

Kepler’s equation is then solved for E in an interative procedure as 

E n = M + e sin E n _ t 

where E 0 = M and M = n t 

The true anomaly 6 at time t is then found from 

Vl - e 2 s in E 


(13) 


(14) 


(15) 


sin 6 


cos Q - 


1 - e cos E 

cos E-e 
1 - e cos E 


> 


(16) 




We now have S - 6 -Q and using the procedure in Appendix B we may find 
t , v the position and velocity at time t . 

n n 

Define a function F 

F = D - S (17) 

where D is the distance from Jupiter which is obtained from r n and the ephem- 
eris of Jupiter. A new At n+ x is then calculated from 

At , F - At F . 

n— 1 n n n-1 


At 


n + 1 


where of course at the beginning 


and 


F -F , 

n n-1 


At o - 0 


( 18 ) 


F„.i = - S - 
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We then return to equation (14) and repeat the procedure until the function F is 
within a predetermined small value. 

The Undisturbed Post Encounter Trajectory 

We have now determined the heliocentric position and velocity r and v of 
the probe when it enters Jupiter’s sphere of influence. The next problem is to 
find the post encounter orbit of the probe after it has passed thru the sphere of 
influence as a function of encounter geometry. References 8 thru 13 describes 
various methods used to obtain the post encounter orbit as a function of en- 
counter conditions. The method described below is that used by the authors of 
references 8 and 9 with some modification to make the method more accurate. 

We have v from above and knowing the injection date of the probe and the time 
to the sphere of influence we can find v p the heliocentric velocity of Jupiter 
from the ephemeris data. The velocity of the probe with respect to Jupiter at 
entrance to the sphere of influence, v J , is then given by 



as can be seen from a diagram of the velocity vectors. 

Within the sphere of influence the probe's trajectory will be hyperbolic 
with respect to the planet (reference 1 — see figure 5) . We assume that we 
may choose the radius of closest approach r of the probe without significantly 
altering the direction or magnitude of v' . This is equivalent to assuming that 
while the direction and magnitude of vj are fixed at the calculated values, we 
are still free to choose the point at which the probe pierces the surface of the 
sphere of influence (point A in figure 5). That this assumption is valid is indi- 
cated by Table 1 (from reference 8) where v' lies along the S° vector (figure 
5). This table was compiled for both retrograde and posigrade flybys of Jupiter 
and for a radius of closest approach r p varying from 100,000 to 10,000,000 km. 
These values for r p span a far larger region than we are interested in. 
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Figure 5-Encounter Orbit Within the Sphere of Influence 
(From Reference 8) 



For a given asymptote and energy at entrance it will be shown that the entire 
trajectory within the sphere of influence, including r p , is determined by the 
location of point A. 

A convenient method for specifying the location of point A (reference 8, 9 
and 11) and the radius of closest approach is to select the magnitude of the miss 
vector ( |b| = B in figure 5) and an angle 'P , measured between B and a fixed 
vector in the planet's orbital plane. This fixed vector T° is one axis of an 
orthogonal coordinate system whose origin is the planet's center and consists 
of the unit vectors S° , T° and R° (References 10, 11, 14 and 15). These 


vectors can be defined in the following manner: 
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Table I 


The Velocity Vector as a Function of Choosing the Entrance 
Point on the Sphere of Influence ( from reference 8) 


Posigrade 
r p x 10 S km 

s° 

IvJ 

(km/sec) 

Right Ascension 
(DEG) 

Declination 

(DEG) 

1 

161.38295 

-.037895218 

13.863505 

5 

161.36287 

-.037826118 

13.873625 

10 

161.34038 

-.037764311 

13.882193 

50 

161.10998 

-.037458739 

13.910651 

100 

160.75975 

-.037179910 

13.922712 

Retrograde 




1 

161.43148 

-.037922481 

13.863560 

5 

161.48549 

-.037896001 

13.873828 

10 

161.53529 

-.037877017 

13.882264 

50 

161.82623 

-.037871506 

13.910774 

100 

162.11483 

-.037958746 

13.923062 


S° is assumed to lie along the incoming asymptote of the hyperbola; 

T° = k° x S° 


where k° is a unit vector perpendicular to the planet’s orbital plane; 

R° = S° x T°. 


The vector T° lies in the planet's orbital plane and is positive in the direction 
of the planet’s motion. We define W to be the angle between T ° and B measured 
in the direction T° to R °. Figure 6 (From ref. 14) shows the geometry of this 
coordinate system and its relationship to the impact plane which is perpendicu- 
lar to the incoming asymptote. 




16 



Figure 6-Planetary Approach Coordinate System and Impact Plane (From Referenced) 


The total effect of the undisturbed encounter is simply to rotate , thru 
an angle y , as illustrated in figure 5, since we must have 



if energy is to be conserved (references 1, 8 and 11). We may easily calculate 
y if B is specified: 


let a be the semi major axis of the hyperbolic orbit, 
straight forward rearrangement of the energy equation we 


s 


2U -S|v; 


' I 2 


Then from a 
have 


a 


(19) 
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where /j. is the planet's gravitational constant. From figure 5 and the proper- 
ties of an hyperbola we have the eccentricity 

e = _J_ (20) 

cos e 


where e is the half angle between the asymptotes and 

B 

sin e = = — cose. 

a e a 

Therefore 


Tan e = — (21) 

a 


now 

7 = 77-26 or 

y = 7T - 2 Tan" 1 (22) 

and the radius of closest approach is 

r = a (1 - e ) . (23) 

We have already pointed out how we may choose the point of entry on the 
sphere of influence without altering the direction or magnitude of vj . However, 
this obviously will change r ' the Jovian centered position vector of the probe. 
We have shown how the parameters B and T define the entrance point on the 
sphere of influence and allow the calculation of the hyperbolic orbit. We can 
use B and w to obtain r x from which we will then be able to easily calculate 
the post encounter orbit at exit from the sphere of influence. We will also use 
v ' and the resulting r J as starting conditions for the calculation of the disturbed 
orbit in the next section. 

Figure 7 shows the relationship of the vectors and angles involved in the 
R°, S°, T° coordinate system. 
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Figure 7-Relationships in the Approach Coordinate 
System (From Reference 8) 


Referring to this figure we see that 

B° = T° cos ¥ + R° sin 'P 

From equations (19) and (20) we find the semi latus rectum 

p = | a (e 2 - 1) | (24) 

and the true anomaly at entrance to the sphere of influence 

6 - cos" 1 — (25) 

Se 

Now from figure 5 we find r ^ , the Jovian centered position vector at entrance 
to the sphere of influence 

r ; = S [cos (0 + e - 77/2) B° + cos (< 9 + e)S°] (26) 


Equation (26) however does not represent the true position vector since it was 
assumed that v \ represents the incoming asymptote 
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S° is therefore not the true asymptote and the miss vector B is not quite 
perpendicular to S°. The position vector r ^ is thus not quite compatible with 
the anomaly 6 . We may make the entire orbit compatible by using r'j to cal- 
culate a new 6 and a new time in the sphere of influence as follows: 
the semi latus rectum is 

2 


P = 


r i* v i 


then 


C3 = 


Vn 


a = - 


and finally 


C3 


-f* 


C3 p 


6 - cos 


-l 


' p - V. 


r l e 


The procedure in Appendix A then yields the time from perijove to the sphere 
of influence which is one half the total time spent in the sphere. The differences 
in total time spent in the sphere of influence between this new trajectory and 
the approximate trajectory is on the order of 10 minutes when the total time is 
approximately 80 days and the magnitude B of the miss vector changes by less 
than one percent. Since B was one of the parameters which define the en- 
counter orbit and therefore is chosen beforehand, any change in B will slightly 
alter the encounter trajectory. But since the change in B is small any changes 
produced by making the orbit consistent with r J are negligible. Table 2 shows 
the differences, resulting from this procedure, obtained in the magnitude of the 
miss vector and the total time in the sphere of influence. 
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Table II 

Effects of Approximating the Encounter Trajectory 



Old B (km) 

New B (km) 

A (km) 

Old Time (days) 

New Time (days) 

A (days) 

180° 

1.25 x 10 6 

1.26 x 10« 

.01 x 10 6 

83.009384 

83.017077 

.0077 

0° 

1.25 x 10 6 

1.26 x 10 6 

.01 x 106 

83.009384 

83.017128 

.0077 

180° 

1.0 x 10 6 

1.0078 x 10 6 

.0078 x 10 6 

81.453790 

81.460214 

.0064 

0° 

1.0 x 10 6 

1.0078 x 10 6 

.0078 x 10 6 

81.453790 

81.460217 

.0064 

180° 

5.0 x 10 6 

5.0389 x 10 6 

.0389 x 10 6 

82.912777 

82.915892 

.0031 

0° 

5.0 x 10 6 

5.0388 x 10« 

.0388 x 10 6 

82.912777 

82.915926 

.0031 


With the now determined Jovian centered position and velocity r [ and v ' 
and the total time t spent in the sphere of influence we are in a position to find 
r ' and V 2 the Jovian centered position and velocity of the probe at exit from 
the sphere (Point A' in figure 5) . To do this one needs an accurate means of 
stepping along the hyperbola. The method used to determine r ’ and v' is 
described fully in Appendix C. 

Using t the time in the sphere of influence we can determine the position 
and velocity of Jupiter from ephemeris information at the time the probe leaves 
the sphere of influence. This position and velocity with r* and determines 
the heliocentric position and velocity of the post encounter trajectory. 












THE PERTURBED JUPITER FLYBY 


As pointed out in the previous section the Jovian centered r ' and v [ 
obtained there are also to be used as starting conditions for the disturbed tra- 
jectory. The disturbed trajectory is treated in the standard manner using 
either Cowell’s or Encke’s method. A review of these methods is given below. 

Cowell’s Equations of Motion - Perturbations in 
Rectangular Coordinates 

Figure 8 shows three bodies, S , P and P ' with masses M , m and m' re- 
spectively. The point 0 is a fixed point in space. We consider P to be the 
body which is in motion about S and is disturbed by body P' . 


P(m) 



0 


Figure 8-Vector Diagram of Cowell’s Equations 


21 



22 


Referring everything to point o we may express the force on S as 
^ r * k 2 Mm — -* - k 2 Mm — — 


M 


dt 2 

The force on P is 
md 2 7 


|r x - r 2 l 


( r 2 - r,) +- 


( r 3 ' r l> 


(27) 


r l - r 3 


dt : 


2 k 2 Him — * v k 2 m m 1 — 

~ = ri — ( r i - r 2> + "~ TT (r 3 - r 2 > 


(28) 


|r 2 — r i! 


I r 3- r 2 ! 


we define 



r ' - r = r 3 - r 2 - p 

By subtracting equation (27) from (28) and making the above substitutions we 
obtain 


d 2 7 


dt 2 


k 2 (M + m) _ , , 

r + k^ 


m 




(29) 


This is Cowell's equation of motion in rectangular coordinates and is now 
referenced to the central body S . If there is more than one disturbing body 
the above equation may be written as: 


d 2 7 
dt 2 


k 2 (M + m) 



1 


m . 




(30) 


Encke's Equations of Motion 

In the above development if the acceleration due to the disturbing bodies is 
of the same or higher order as that due to the central force field then Cowell's 
procedure is a good one to use (reference 1). If, however, the disturbing ac- 
celeration is small then Cowell's method becomes inefficient. If in this case 
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only the differential accelerations are integrated then efficiency and accuracy 
may be retained. This procedure is the well-known Encke's method. 

Define 

7 = 7 0 + S r 


where 7 0 is the two body position vector and Sr is the differential correction. 
The equation of motion is then 


d 2 7 _ k 2 Mr 
dT 1 ~ 


+ F 


(31) 


where F is the disturbing acceleration. 
For the two body case we have 


d 2 r Q k 2 M r 0 


dt 2 r * 


(32) 


Upon subtracting (32) from (31) we obtain 


d 2 S7 k 2 M8 7 , 2 .. / 1 1 


+ k 2 M / ] r + F 


dt 


A r3 


Denote 


k 2 M 


= h 


(33) 


(34) 


and 


— = 1 + 2q 


Observe that since r 2 and r 2 are nearly equal, 2q will be small. 


1 + 2q = 


7-7 ( 7 o + + s7 > 


_ r 0 2 + 27 0 • 87 + S7 • 87 


Now 
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Therefore 


q = 


(r Q + 1/2 Sr) • Sr 


(35) 


and 


i-=(l + 2q)- 3 ' 2 

r 3 


Then equation (33) becomes upon substituting (34) and (36) 
d 2 SI 


Observe that 


then 


dt : 


- = - hSr + h [1 - (1 + 2q)" 3/2 3 r + F 


_3/2 /, 15 \ 

1 - (1 + 2q) = q \ 3 q +••*] = 


fq 


d 2 Sr 


dt ' 


— hS"? + hf q7 + F 


(36) 


(37) 


(38) 


Equation (38) is the form of Encke’s equations used for computation. However, 
Encke’s equations do have a disadvantage that causes computation time to in- 
crease. When the magnitude of Sir in equation (38) becomes large it is neces- 
sary to change the two body vector r 0 by defining r 0 = r if accuracy is to be 
maintained. This process, known as rectification, requires the numerical in- 
tegrator to be restarted with a resulting increase in time. It should be pointed 
out that the disturbing acceleration F is the same as in equation (30) above; 
i.e., 


V — i / r ' - r r' 

= > k 2 m! _L_ 

“ 1 V p 3 r '. 3 

i \ r i i i 


( 39 ) 


and that the only change has been in the central body term of equation (30) . 
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Perturbations Due to the Sun 

In the above equation (39) for the disturbing acceleration we will include 
only the four great moons of Jupiter. Since the distance of the sun from Jupiter 
is so much larger than that of the moons, the term in parentheses in equation 
(39) for the sun will lose accuracy because of the fact that both members are 
almost equal. We may eliminate this problem in a manner similar to Encke's 
procedure by treating the sun as a special case. Reference 16 describes a 
method for solving this problem and the procedure is repeated below. 


JUPITER 



PROBE 


Figure 9 shows the geometry of the sun, Jupiter and Probe during the 
encounter phase. By a slight change in notation we may write equation (29) for 
the case where Jupiter is the central body and the sun is the disturbing body: 

dt 2 r 3 \ | (r ) + r| 3 |(r) | 3 / 

or 

+ ? (4( 

dt 2 r 3 \|(r) + r 1 3 |(r) | 3 / 


where r and (r) are defined as in figure 9 and the notational change is obvious 
by comparing figures 8 and 9. 
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Now 


1(7) + r | 2 = |(7)| 2 + 2(7) • 7 + |7 


(41) 


let 


then 


q = 


-1- r 2 + (7) • 7 


I (7) I ' 


2 | (7) | 2 q = r 2 + 2 (7) • 7 
upon substituting this into equation (41) we obtain 

|(7) +71 2 = 1(7) I 2 + 2 I (7) I 2 q 
= I (7) | 2 (1 + 2q) 


then 


let 


|(7) + 7| " 3 


1 .. ~ v -3/2 

(1 + 2q) 


1(7)1 


h = 


k 2 m2| 


y 


(h) = 


k 2 M 

1(7)1 




Substituting equations (43) and (44) into (40) yields 

— = -h7-k 2 M J((7)+ 7)_ L_ (l + 2q)" 3/2 - 


dt 2 


1(7)1 


1(7)1 


(42) 


(43) 


(44) 


= - hr - (h) [(7) +7] (1 + 2q) “ 3/2 + (h) (7) 


— = [_h -(h) (1 + 2q)" 3/2 ] 7 + (h) (7) [-(1 + 2q) _3/2 + 1] 
dt 2 


( 45 ) 
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now making the same series expansion as before 

1 -(l + 2q)‘ 3/2 = q^3 -^-q +^.q 2 + • • 

= fq 

the equations of motion become 

= {-h- (h) [1 - fq] > r + (h) (r) fq 

dt 2 

or 

— = {(h) fq - (h) - h} r + (h) fq (?) (46) 

dt 2 

where the term due solely to the central body is hr and the rest is due to the 
presence of the sun. 


Combining the Equations of Motion 

The equations of motion are now combined to include the perturbations due 
to the sun and those due to other bodies. Rewrite equation (46) in the following 
form 

h 2 r 

= - hr + 4 > (47) 

dt 2 

where 4 > is the acceleration due to the presence of the sun 

4 > = {(h)fq - (h)} 7 + (h) fq (7) (48) 

and 

k 2 mj, 

h = — 

r 3 

where r is the distance from Jupiter to the probe. 
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To include other bodies we must add to equation (47) the acceleration expressed 
as F in equation (39) 

— = - hr +4> + F (49) 

dt 2 

where 



and r is the vector from Jupiter to the probe and r ' is from Jupiter to the dis- 
turbing body. 

Equation (49) represents Cowell's equation of motion. When Encke’s pro- 
cedure is wanted the following is used: 


d 2 S: 


k 2 M. k 2 m^ 
2 Sr + — 


fq r + F + <p. 


(50) 


dt 2 r o r o 

Care must be taken in interpreting equations (48), (49) and (50) since f and q 
have different meanings depending on whether or not Encke’s method is used. 


Integration of the Equations of Motion 


The above equations of motion either (49) or (50) are integrated over the 
time spent in the sphere of influence by means of a numerical integration 
technique (reference 17). This technique is a fourth order predictor- corrector 
with automatic adjustment of the time step of integration. 

The integration is started at the time of the probe’s entrance to the sphere 
of influence and is terminated at exit from the sphere. This termination is 
accomplished by checking the probe’s distance from Jupiter until this distance 
becomes equal to the radius of the sphere. The integration is then halted and 
the resulting trajectory is determined. 
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In the actual computer runs made, the use of the above Encke's equation 
(50) more than tripled the computer time when the moons of Jupiter were in- 
cluded as perturbing bodies. This increase in time was due in part to the 
greater number of calculations required and to the number of rectifications that 
occurred because of the inclusion of the moons. For this reason Cowell's 
method was used exclusively for the results included in this report. Encke's 
method was used as a check in only a few cases to insure that numerical ac- 
curacy was being retained; 



THE GALILEAN SATELLITES OF JUPITER 


On the night of January 7, 1610 Galileo Galilei turned his crude telescope 
toward Jupiter and noticed three small "stars” in the vicinity of and in a straight 
line with the planet. Repeated observations showed that these "stars" moved 
with respect to Jupiter and on January 13 he saw four "stars". Further observa- 
tions proved to him that he had indeed discovered satellites of Jupiter (reference 
18). 

The German astronomer Marius later independently discovered the Galilean 
moons and gave them their names (reference 19) as follows: 

I Io 

II Europa 

DI Ganymede 

IV Callisto 

Table III is a compilation of the physical data of these satellites and is 
taken in part from reference 20. 

When we consider that the planet Mercury has on the same scale a diam- 
eter of 0.38, compared to 0.394 for Ganymede, we see that these moons are by 
no means small. 

These four satellites have been found to be subject to very strong mutual 
perturbations which cause librations in their orbits about Jupiter. Laplace 
first discovered a resonance relationship between satellites I, n and III which 
shows that their mean motions are approximately in the proportion of 4:2:1 or 
in terms of mean longitudes L 1 -3L 2 +2L 3 =180 degrees (reference 21). If one 
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Table III 

The Galilean Satellites: Astronomical and Physical Data 



h 

Europa 

Ganymede 

Callisto 

Mass (Earth = 1) 

0.0121 

0.0079 

0.0261 

0.0160 

Mean Diameter (Earth = 1) 

0.255 

0.226 

0.394 

0.350 

Mean Density (H 2 0 = 1) 

4.03 

3.78 

2.35 

2.06 

Surface Gravity (Earth = 1) 

0.19 

0.16 

0.17 

0.13 

Velocity of Escape (km/sec) 

2.5 

2.1 

2.9 

2.4 

Period of Axial Rotation (Days) 

1.77 

3.55 

7.15 

16.69 

Mean Distance from Planet (km) 

421,400 

670,500 

1,069,500 

1,881,200 

Sidereal Period (Days) 

1.769 

3.551 

7.155 

16.689 

Mean Orbital Velocity (km/sec) 

17.37 

13.77 

10.90 

8.22 

Radius of Sphere of Influence (km) 

7,200 

9,600 

24,800 

36,100 


of the satellites is disturbed, the motions of the other two adjust so that this 
relationship is satisfied. Therefore the three satellites cannot be in conjunction 
or opposition at the same time (reference 22). 

Little can be found in the literature on an accurate theory of the motions of 
the four moons but accurate positions are obviously necessary if a realistic 
investigation into their effects on the flyby probe is to be made. The theory of 
the moons is quite complex and interest in the theory has only recently begun to 
bring about useful investigations. 

Sampson in 1910 (references 23 and 24) and De Sitter in 1931 (reference 
25) have developed detailed theories to describe the motions of the Galilean 
satellites. Sampson's work is perhaps the best known and includes a rather 
extensive set of tables which are used to determine the position of the moons 
and to predict such physical phenomena as occulations and eclipses. This work 
which has become known as Sampson's Tables represents a monumental and 
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complex task. Although difficult to use these tables give the times of phenomena 
to 0.000001 days and positions to 0°. 000001 (reference 26). Andoyer (reference 
27) somewhat simplified the procedure by using only the main terms of Samp- 
son's tables and even though this is less accurate, Andoyer's procedure has 
been adopted by the Nautical Almanac Office in their publication of the pheno- 
mena of the moons. 

More recently Ferraz Mello (reference 28) has studied the planar motion of 
the Galilean satellites in rectangular coordinates. He shows that the higher 
order harmonics of Jupiter's gravitational potential and the relativity correc- 
tions are negligable. Marsden (reference 29) has studied the motion of the four 
moons using Von Ziepel's method to eliminate the short period terms and to 
reduce the number of degrees of freedom. 

The method used in this report to obtain the positions is a more or less 
simplified version of Sampson's method. The model used is one of a slowly 
moving plane common to all four satellites with some plane perturbations. This 
model leads to uncertainties in position of 

300; 6000; 6000; 25000 kilometers 

respectively for satellites I, II, III and IV when compared with Sampson's Tables 
(reference 30). This method was received from Dr. Jean L. Sagnier in Paris, 
France who has spent much time investigating the theory of the moons of 
Jupiter. 

A complete description of the information received from Dr. Sagnier ap- 
pears in Appendix D and a complete and extensive discussion of Sampson's 
theory of the moons can be found in reference 23. The Tables in Appendix D 
are used to determine the positions of the moons which are in turn necessary to 
evaluate the disturbing acceleration in equation (39) above. This method of 
obtaining the positions of the moons has been programmed to generate an 
ephemeris of positions for each of the four moons. Such ephemerides do not 
seem to be easily available elsewhere. 
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A comparison was made with examples worked out by Sampson in reference 
20. Table IV shows the results of this comparison for satellites I and in on two 
different dates. The longitude is measured from the mean equinox of 1900 along 
the mean ecliptic to the ascending node of Jupiter’s orbital plane and there- 
after along this plane. The magnitude of the radius vector to the satellite is 
the projection of the radius vector into Jupiter’s equatorial plane; and the 
column labeled "tan latitude" is the tangent of the elevation of the radius vector 
to Jupiter's orbital plane. 


Table IV 

Comparison of Sagnier’s Method and Sampson's Tables (Jovicentered) 


Satellite I June 2.156 1909 
simplified method 
Sampson's tables 
I difference I 

Satellite in April 3.841 1910 
simplified method 
Sampson's tables 
I difference I 


Longitude 
336°. 07021 
336°. 06863 
0°. 00158 


Radii (A. U.) 


0.0000281304 
0.0000281296 
8 x IQ" 10 


Tan Latitude 
.018832 
.0191815 
.0003 


9° 


9° 


0 ° 


.92885 

.92550 

.00335 


0.00714796 
0.00714794 
2 x 10- 8 


.045659 

.045686 

.000027 
















ENCOUNTER GEOMETRIES AND MISSIONS USED IN THE STUDY 


Twenty different trajectories were studied to determine the effects of the 
Galilean moons. The launch date for injection into the transfer orbit from the 
Earth to Jupiter was March 8, 1972 which is one of the dates under considera- 
tion for a Galactic Jupiter Probe currently being planned by NASA (reference 
31). Flight times of 450, 500, 550 and 600 days from the Earth were used along 
with several combinations of the encounter parameters B and V . The first 
twelve cases included flight times of 450, 500 and 550 days. The magnitude of 
the miss vector was held constant at 1.25 x 10 6 km and the angle ^ was varied 
from 0 to 360 degrees in steps of 90 degrees. The remaining 8 cases were 
obtained using a constant value of 180° for and values of 1.0 x io 6 and 1.5 x 
10 6 km for B . A flight time of 600 days was also included in these cases. 

These 8 cases represent the type of encounter geometry that is particularly 
useful in obtaining hyperbolic heliocentric trajectories which are capable of 
achieving far distances in relatively short times. 

Some of the first twelve trajectories, particularly those where 'V = 0° , are 
also useful for solar probes since they result in eccentricities of almost 1. 

These trajectories can in some cases impact the sun which cannot be done with 
direct flights from the Earth because of the energy required (references 8 and 
11). However for these trajectories it is desirable that at exit from Jupiter's 
sphere of influence the probe be heading toward the sun, i.e. that the heliocentric 
true anomaly be greater than 180°. Otherwise the time required for the probe 
to reach the sun is increased prohibitively. 
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Values of 90° and 270° for 'I* are very useful in obtaining large distances 
both above and below the ecliptic plane and a little reflection will reveal that 
the two values ought to give symmetrical results (reference 8) . Indeed it is 
possible to obtain a distance of 30 a.u. out of the ecliptic as did the case for 
450 days where V = 270°. However the time required to reach such distances 
will also be considerable and may rule out the possibility of such a flight being 
undertaken. 

Other types of missions may also be studied if desired. These could in- 
clude obtaining high inclinations of the orbit plane with the ecliptic, maximum 
velocity normal to the ecliptic and many other interesting orbits. Most of these 
would be difficult to attain in terms of the energy requirements without using 
the gravitational attraction of a large planet. This report includes only three 
types of mission objectives: the time required to reach 10 a.u. from the sun, 
the maximum distance attained out of the ecliptic, and close approaches to or 
impact with the sun. 


I 



RESULTS AND CONCLUSIONS 


The effects of the Galilean moons of Jupiter on the Jovian centered orbit 
are presented in Table V. These results are presented as differences resulting 
from using only the sun as a disturbing body and then using the sun and moons 
as disturbing bodies. The first three columns identify the case with the param- 
eters B in kilometers, W in degrees and transfer time in days. The next six 
columns are the differences in the Keplerian elements as follows: the semi 

major axis (a), eccentricity (e), inclination (i), true anomaly ( 6 ), right ascen- 
sion of the ascending node (D), and arguement of perojove (^). The next column 
is the difference in magnitude of the velocity at exit from Jupiter's sphere of in- 
fluence. As described in a previous section the integration of the equations of motion 
is terminated when the probe reaches a distance of 48 x 10 6 kilometers from 
Jupiter which is the radius of the sphere of influence. Therefore, the magnitude 
of the radius vector is always the same at exit and differences can occur only 
in tangential position and in the velocity vector. The final four columns show 
the closest approach to the Galilean moons. The closest approaches to all four 
moons always occur within a span of 1.5 days and always near perijove of the 
approach trajectory where the probe is moving very fast. 

A quick glance at the A a column shows that the only large differences 

occur in cases 11 and 19. The semi major axis is computed from 

r 

a = 
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and since r is constant, the change in a is dependent only on v. Case 11 
passes Europa by 6100 km and is the closest approach of all the cases. A close 
approach to one or more of the moons will of course perturb the trajectory. A 
study was made for case 11 in which one moon at a time was included into the 
equations of motion. Europa by itself caused about half of the perturbation and 
inclusion of Ganymede accounted for almost all of the remaining perturbation. 
However, Ganymede by itself accounted for little of the total perturbation. Hand 
calculations independent of the program substantiated the results obtained when 
a single moon was included. This study shows that a complex dynamical inter- 
action between the bodies is occurring. This complex interaction cannot be 
predicted analytically in an easy way and therefore the perturbation is difficult 
to prove conclusively. 

The change occurring in case 19 is even more difficult to explain. This 
trajectory does come closer to Ganymede than any other trajectory but is well 
outside Ganymede’s sphere of influence (see table III). This case represents a 
600 day transfer trajectory and as such is a relatively low energy Jupiter 
centered trajectory. Because of its low energy it is more sensitive to perturb- 
ing forces. As in case 11 approximate methods tended only to confirm the re- 
sults obtained using one moon at a time. Again the explanation lies in the 
dynamical interaction between the bodies. 

The other 18 cases do not exhibit any abnormally large changes in any of 
the elements. 

It therefore appears from Table V that the moons of Jupiter do not repre- 
sent any major perturbative force as far as the Jovian centered trajectory is 
concerned in a majority of the cases. This is consistent with what was originally 
thought when the problem was first undertaken. However, cases 11 and 19 
indicate that under the right dynamical conditions, significant perturbations may 
result that will have to be considered. These perturbations cannot be predicted 
without including the moons. 
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Table VI shows the resulting effects of the moons on the heliocentric tra- 
jectories at post encounter. The columns are the same as for Table V except a 
column has been added for the difference in the heliocentric radial distances. 

Large differences occur in the semi major axis for the heliocentric orbits, 
which by virtue of the above equation are attributed to changes in the position 
and velocity magnitudes. Although these are large variations, it must be pointed 
out that the semi major axis is itself large. For example Case 12 shows a 
change in a of .12 x io 9 km, where a equals 3 x 10 9 km, representing a per- 
turbation of about 4%. This is the largest percentage of all 20 cases. 

Cases 11 and 19 in the heliocentric orbits do not exhibit the large perturba- 
tions that occurred in the Jovicentric orbits. The variations in the semi major 
axis in the heliocentric cases are much less than 1%. Indeed a glance at Table VI 
shows that none of the elements change very much and some exhibit smaller 
changes than the other 18 cases. 

Table VTI shows the differences resulting in the parameters that are of 
interest for a deep space post Jupiter encounter mission. The first three 
columns are the same as the previous tables. Columns 4 and 5 are found by 
taking the square root of the sum of the squares of the differences in com- 
ponents of the Jovicentric position and velocity vectors respectively. These 
numbers represent the straight line distance and velocity between the two types 
of perturbation rims. Column 6 is the difference in the time spent in Jupiter's 
sphere of influence. Columns 7, 8 and 9 are the respective differences in time 
to reach 10 a.u. from the sun; distances attained out of the ecliptic; and closest 
approach to the sun where these missions are applicable. Again the last 4 
columns are a repetition of the close approaches to the moons. 

Figure 10 (from reference 9) is a contour map of constant times to 10 A.U. 
This map was done for a 500 day transfer time and is in the T and R plane 
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which is defined in a previous section. This plane is commonly called the impact 
plane and is perpendicular to the incoming asymptote. The angle 4* is measured 
around the perimeter of the figure and the magnitude of the miss vector (B) is 
measured by the concentric circles. A flight time of 600 days from post en- 
counter to 10 A.U. can be achieved by any combination of the parameters B 
and 4* as obtained from the 600 day contour line. Anything inside this line rep- 
resents a flight time of less than 600 days. The launch date for this figure was 
March 10, 1972 which is two days later than the launch date for Case 7 in 
Table V, but this difference will not alter Figure 10 by any appreciable amount. 
Figure 10 was drawn using only a two body determination as outlined in a pre- 
vious section. Case 7 in Table VII took 558.990 days to reach 10 A. U. when the 
sun and moons were included as perturbing bodies. The undisturbed time in 
Figure 10 was 556.528 days. From Table VII the time difference caused by the 
moons of Jupiter was 0.12 days. Cases 3 and 4 in Table VI may also be located 
on this contour map for B equal to 1.0 x 10 6 and 1.5 x 10 6 km respectively. 
Again it is seen that the influence of the moons is negligible. Also the dif- 
ference between the perturbed and two body determination used in Figure 10 
was less than four days so even the sun makes only a small difference as a 
disturbing body. 

As far as the mission is concerned it is seen that the effects of the moons 
are very small indeed and can be neglected. However as far as determining 
the precise point of exit from the sphere of influence, columns 4 and 5 of 
Table VII clearly indicate that the moons do have an effect. These effects could 
be of critical importance depending on the particular post encounter mission 
being contemplated. This along with the changes in the heliocentric semi-major 
axis is sufficient to conclude that for purposes of orbit determination the 
Galilean moons must be included in any realistic study of a Jupiter Flyby. 


Effects of the Galilean Moons on the Jovicentric Orbital Elements 
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Effects of the Galilean Moons on Deep Space Mission Parameters 
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SUMMARY 


The effects of the Galilean moons of Jupiter on a probe passing close to 
Jupiter were studied by using two different perturbative schemes. The first 
included only the sun as a perturbing body and the second included the sun and 
the four moons. Variations obtained from the two models are then the effects 
of the presence of the moons. An integrated trajectory was used only within 
Jupiter's sphere of influence with an approximate two body determination of the 
transfer trajectory between the Earth and Jupiter. 

The numerical integration method was a predictor- corrector scheme with 
automatic adjustment of the step size and automatic error control. Due to con- 
siderable savings in machine time all cases were studied by integrating Cowell's 
equations of motion rather than Encke's equations. 

While the total variations due to the Galilean moons in the Jovian centered 
orbital elements of the probe during the Jupiter flyby are fairly small, the 
actual point of exit on the sphere of influence is changed by an appreciable 
amount. Therefore as concerning gross trajectory determination one can es- 
sentially neglect these moons but when accurate position and velocity vectors 
are required, the perturbative effects of the moons most certainly should be in- 
cluded. 

It is shown that the post encounter mission objectives such as time to reach 
10 a.u. from the sun are also little changed due to the presence of the four larger 
moons of Jupiter. However the semi-major axis of the resulting heliocentric 
orbits do undergo large variations. The importance of the perturbations will 
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depend on the exact nature of the post encounter missions that are to be studied 
such as a close approach to a second outer planet. Thus for accurate trajectory 
determination as well as accurate positions and velocities it is advisable to in- 
clude the Galilean moons in the equations of motion. 


APPENDIX A 


TO COMPUTE TIME FROM PERIAPSIS TO A GIVEN TRUE ANOMALY 

This problem is easily solved using two body equations and the solution of 

Kepler’s equation yields the desired result immediately. 

Let: a be the semi-major axis 
e the eccentricity 

6 the true anomaly 

n the mean motion 

p the semi latus rectum 

t the time from periapsis 

m the mean anomaly 

E the eccentric anomaly 

r magnitude of radius vector 

The relationship between 6 and E can be seen from figure A-l. 

Now X = E, - ae = a (cos E - e) = r cos d (A-l) 

Y = r] = bsinE=r sin 0 
b = a Vl - e 2 

also 

a (1 - e 2 ) 
r = — ^ 

1 + e cos 6 

from (A-l) and (A-2) we obtain 

r = a (1 - e cos E) 


(A-2) 


(A-3) 


(A-4) 
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from (A-l), (A-2) , and (A-4), we obtain 


cos 6 - 


cos E - e 
1 - e cos E 


sin 


e = 


Vl - e 2 s in E 
1 - e cos E 


and this allows us to finally obtain 




(A-5) 


The mean motion n is defined to be 



where /z is the gravitational constant of the central body. Using (A-5) to obtain 
E , Kepler’s equation states that for an ellipse 


M = E - e sin, E = nt 


Therefore 

t=“ 

n 

For an hyperbola, Kepler’s equation becomes 


(A-6) 

(A-7) 


M = e sinh E - E 


(A- 8) 


and again equation (A-7) holds. 

For a parabola, a slightly different treatment is necessary. Equation (A-3) may 


be written as follows: 
let q = a(l-e) 


r _ a(l -e) (1 + e) 
1 + e cos 0 


then upon setting e = 1 


Define 


then 


2q 

r = 

1 + cos 6 

D = V2q tan| 
r =1 (2q +D 2 ) 


(A- 9) 


(A- 10) 


I 
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The area using equation (A-10) is 


This yields 


where 



n-iJT 


(A- 11) 


equation (A-ll) may be found developed fully in reference 3. Equations (A-5), 
(A- 6), (A- 7), and (A-ll) then allows one to compute the time for any type of orbit. 



Figure A-l — Relationsh ip Between True Anomaly and 
Eccentric Anomaly 


APPENDIX B 


TO UPDATE THE CARTESIAN COORDINATES OF THE PROBE 
USING AN INCREMENTAL TRUE ANOMALY 

The problem to be solved here is that of determining the position and 
velocity of a probe which is separated from a known position and velocity by an 
increment in the true anomaly. This treatment is essentially that found in 
reference 7. The following notation will be used: 

(j. gravitational constant of the central body 
r o initial position 
v 0 initial velocity 
H = r 0 x v Q momentum vector 

b = I H | 

S = 6 -d Q the incremental true anomaly 
r final position 
v final velocity 

Since we are considering only two body planer motion we may write the 
following relationships: 

7 = f 7 0 + gv I 

f (B_1 

7 = f 7 0 + g v 0 J 

The quantities we wish to find are f , g, f, and g. 
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Comparing this with equation (B-l), we see that 


= rr ( c 


£ XI 0 0 

f = — cos S « — si 


0 

r r , 


in 


g = 1 sin S 

6 h 




(B-4) 


J 


However, the magnitude r in the above equation is still unknown. To solve for 
r we make use of the radius equation 

r = a < 1 ~ e2 ) = 


1 + e cos 6 1 + e cos 6 


and 


Thus 


also 


Therefore 


a(l - e 2 ) = 




r = 


/x(l + e cos 6 ) 


9 = s + e, 


cos 9 = cos 8 cos 6 0 - sin 8 sin 9 Q 


Using equations (B-5) and (B-6), we immediately obtain 

h 2 - flT 0 


cos 6 0 = • 


n e 


Note that from the radius equation we may write 


* d ^ o v o 

r ° = dt (r » ' r o> = — 


epd sin 6 
(1 + e cos & 0 ) 2 


and by using 


(B-5) 


(B-6) 


(B-7) 


h = r o 
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we obtain 


sin 6 Q = ■ 


h r, 


Me 


from above 


r n • v o 


r o = — f 


which together with equation (B-8) gives 


sin 8 . = 


hr o • v o 

r n M e 


using equations (B-5), (B-6), (B-7), and (B-10) we obtain 


hf 

M 


1 + e 


cos 8 (h 2 - M r 0 ) 


- sin 8 sin 6>, 


M r n e 


(B-8) 


(B-9) 


(B-10) 


or 


M r 0 + (h 2 - M r 0 ) cos 8 - h (Jo ‘ v 0 ) sin 8 


(B-ll) 


Equation (B-ll) allows the computation of f and g in equation (B-4) above. 
In an entirely similar manner we may also obtain f and g: 


f = 


M 


r 0 h 2 


[(r Q • v 0 ) (1 - cos 8) -h sin 8] 


• ^ ** 0 

g = 1 r— (1 - COS 8 ) 


(B-12) 


h S 

We may gain as a dividend the time difference corresponding to 8 in the 
following manner: _ 


a = 


2 - 


r o v ' 


_ h 2 

P~ — - a (l - e 2 ) 
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which yields 


e 2 = 1 - — 
a 

From equations (B-7) and B-10) for 6 we obtain 

h 2 - fXT 

then 


0 = 0 O + S 

and the method of Appendix A will yield the time difference between 0 0 and 6 . 



APPENDIX C 


TO FIND r AND v OF A PROBE ON AN HYPERBOLIC 
ORBIT SEPARATED FROM r 0 AND v Q BY TIME At 

The following method is used to step along an hyperbolic orbit knowing an 
initial position and velocity r 0 and v 0 and a time difference At . 

Reference to figure C-l shows the hyperbola and the geometry of the prob- 
lem to be solved. The axis T, and i are unit vectors, origin at the focus and 

$ v 

such that they form a right handed system. 



Figure C-l— The Hyperbolic Orbit Geometry 
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From figure C-l it is obvious that 

t — > — > 

r = r cos 6 i^ + r sin 6 i ^ (C-l) 

and the well known relationships for a hyperbola are: 


r = a (e cosh F - 1 ) 

(C— 2) 

a e - cosh F 
cos 0 ~ 

e cosh F-l 

(C-3) 

a Ve 2 - 1 sinh F 

sin 9 - 

e cosh F-l 

(C-4) 


where F is the "hyperbolic" eccentric anomaly. Upon substituting (C-2), (C-3), 
and (C-4) into (C-l) , we obtain 


7 = a(e - cosh F) + Yap sinhF i (C-5) 

where p is the semi latus rectum. Differentiating equation (C-5) and noting 
that from Kepler's equation 

dF _ IT 1 

dt f a 3 e cosh H - 1 


we obtain 

— * 7 //,a — * ^ PP — * 

v = £__sinh F + cosh F i ^ (C-6) 

Equations (C-5) and (C-6) allow the computation of position and velocity 
anywhere along the hyperbola if the parameter F is known since a, e, and p 
may be easily calculated from r 0 and v 0 . We will use r 0 and v Q to calculate F 

and to determine the axis i f and i „ . 

s v 

From the Laplace Integral of two body motion we have 


£ 1 - 

E =- v 0 ,c 0 - r o 

where E is a vector in the direction of perihelion and 


c o = r o x v o 


(C-7) 
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Obviously 


also 


*S = 


1 v 


From the energy equation we have 


V E 

c 0 l I2| 


/xr. 


a = 


then 


where 


2 P ~ r o v o 
P = a | (e 2 - 1) | 

e = III 


The mean motion is then 



and the true anomaly at the initial time is 



The change in the "hyperbolic" mean anomaly is 

AM = nAt 


(C-8) 


(C-9) 


(C-10) 


The value for F 0 the eccentric anomaly at the initial point is obtained from 

Ve 2 - 1 sin 6 n 


sinh F 0 = 


e cos 0 0 + 1 


Then the value of the mean anomaly (M) at the point r is obtained by adding 
equation (C-10) to Kepler's equation evaluated at F 0 


M = e sinh F fl - F 0 + AM 


or 


M = e sinhF - F 


(C-ll) 
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where F is the unknown eccentric anomaly needed in equations (C-5) and (C-6). 
Equation (C-ll) can be solved for F by an interative procedure as follows: 
we form the difference function 


A =M-esinhF +F 

n n n 


(C-12) 


and as a first guess we choose F x = 0 . Using the Newton interative formula we 
have 


A 

F , = F - 

n + 1 n 

n 


or 


F 


n + 


= F 


(M - e sinh F n + F n ) 
1 - e cosh F 

n 


(C-13) 


when A n becomes small enough we have the required value for F and r and v 
may now be evaluated. 



APPENDIX D 

Received from Dr. J. L. Sagnier, Paris, France 
THE POSITIONS OF JUPITER'S GREAT SATELLITES 
The following tables include: 

• The initial values and daily motions of the arguments occurring in the 
subsequent series, expressed in radians when taking into account the 
factor 2tt . 

• The series defining the mean positions of Jupiter's equatorial plane and of 
the orbital planes of the satellites. 

• The equatorial true longitudes of the satellites, in radians from 
Sampson's x-axis up to the ascending node of Jupiter’s equator on 
Sampson's (x , y) - plane, and thereafter along Jupiter's equator. 

• The projected values on Jupiter's equator of the radii vectors, expressed 
in A. U. 

• The reductions to Sampson's ( x,y) - plane, to be added to the equatorial 
true longitudes, when the longitudes in Sampson's system are needed. 

• The tangents of the latitudes on Sampson's (x,y) - plane. 

The unit of time is always the ephemeris day, initial epochs are explicitly 
given when necessary. The unit of angle is the radian, except for some phase 
angles in the arguments expressed in degrees. The unit of length is the astro- 
nomical unit. 


58 



59 


A low accuracy determination of the positions of Jupiter’s satellites may be 
conveniently obtained by the introduction of the most significant terms of the 
radii vectors and of the equatorial longitudes in the mean orbital planes, whether 
approximated by the Jovian equator or not. But the use of the series reducing 
longitudes to Sampson's system and of those giving the latitudes in the same, 
is necessary for a full accuracy computation. 


Arguments of Series 


Arg* 

Value for C^TE, 31/12/1959 = 2 tt x 

Daily variation = 

2 77 X 


0.9406 

5442 

8812 

5000 

+ 0.5652 

4720 

1208 

3333 


0.0427 

5197 

0477 

7778 

+ 0.2815 

9656 

0200 

0000 


0.3438 

0071 

0875 

0000 

+ 0.1397 

7123 

9694 

4444 


0.7025 

5215 

8347 

2222 

+ 0.0599 

1974 

8972 

2222 

n i 

0.3517 

7383 

1250 

oooo 

+ 0.0004 

3870 

9722 

2222 

n 2 

0.4138 

9361 

2777 

7778 

+ 0.0001 

3087 

6666 

6667 

77 3 

0.3694 

6197 

9305 

5556 

+ 0.0000 

1930 

9166 

6667 


0.9023 

3018 

4027 

7778 

+ 0.0000 

0527 

0833 

3333 

CO 

1 

0.9339 

4898 

5138 

8889 

- 0.0003 

7230 

6944 

4444 


0.8166 

4930 

4027 

7778 

- 0.0000 

9083 

1388 

8889 

^3 

0.4661 

1489 

5555 

5556 

- 0.0000 

1938 

2222 

2222 

"4 

0.9135 

1122 

8055 

5556 

- 0.0000 

0487 

6111 

1111 


0.8801 

0470 

7361 

1111 

+ 0.0000 

0009 

9722 

2222 

G 

0.6840 

6586 

4444 

4444 

+ 0.0002 

3080 

8888 

8889 

G 

0.5249 

4424 

2500 

0000 

+ 0.0000 

9294 

3888 

8889 

n 

0.0392 

7971 

0277 

7778 

+ 0.0000 

0010 

6111 

1111 


For accurate computations, argument G should be supplemented in each of 
its occurrences by the long period correction: 

2 rr [0. 9500 0000 X 10 - 4 sin 2 tt( 0. 3696 5917 -0.4492 1389 x 10" 4 [t -1960. 0L ) 

v days 7 

+ 0.9208 3333 x 10" 3 s in 2 tt ( 0 . 436 1 4306 +0.2927 7778 x 10" 5 [t -1960. 0] . )] 

v days' 
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Positions of Equatorial and Orbital Planes 

f , 77 .: In Sampson's coordinate system, the first two components of a 
vector of modulus 1, perpendicular to the mean orbital plane of the i th satellite 
(for i = 1, 2, 3, 4) or Jupiter's equatorial plane (for i = 0). 


x 10 6 = 2A sin a, 77 . x 10 6 = X(-A) cos a 
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Equatorial True Longitudes 

For the i th satellite, we have: 

v. = ^ i +2ttx 10"" 8 ‘ I A sin a 

with the following values for A and expressions for a : 



- 

1 097 

1 

to 

+ 39 

M x - M 3 

- 431 

G 

+130 978 

2 V 2 * 2 

+ 6 

“4^ 3 

+ 97 

2G 

+ 

439 

- 3^ 2 

+ 1478 

*1" 77 1 

- 206 

'P-^i 


431 

4 ^- 4 ^ 

+ 1 183 

^1- ^3 

+ 56 

^ -co 2 

+ 

33 

S't 1 -5^2 

+ 528 

^1" 7T 4 

+ 214 

4>-co 2 

+ 

14 

6-1 ^ 2 

- 2 000 

^2 “ 2^3 + 77 1 

- 97 

'P - 

+ 

6 

-U 2 

+ 1 094 

^2- 2 ^3 +77 2 

+ 267 

"2 - 

- 

439 

V <3 

+ 2 517 

V 2 V ff 3 



+ 

369 

2-e,-2X s 

+ 983 

t 2 -2t 3 +7T 4 



SECOND SATELLITE 

- 

3 283 

V ^3 

+ 2 625 

l 2~ 77 2 

+ 94 

2-£ 3 -2G-: 

+297 267 


+ 9 794 

^2“ 77 3 

-2 908 

G 

+ 

1 333 

3^2 - 3^3 

+ 4 119 

^2-^4 j 

- 325 

2G 

+ 

1 503 

4^2 -4^3 

+ 3 169 

^- 2 ^3 + 7r i 

- 236 

5G' - 2G + ^ 

+ 

106 

5^2 -5^3 

+ 2 125 

^2- 2 ^3 + 77 2 

+ 128 

77 3 " 77 

+ 

150 

6^2 -6^3 

-12 375 

' e 2- 2 ^3 + 77 3 

-3 514 

0—0), 


19 

n 2 -n. 

- 4 994 

l 2 - 2 ^3 + ^4 

- 544 

\p - go. 

- 

181 

*2- ^ 

+ 86 

3 ^2 _ 2 ^3 " 77 3' 

+ 256 

0 - 0) ( 

+ 

125 

2^2 “ 2^4 

75 

3^2- 4/t 3 + 7r a 

+ 92 

o 2 - 

+ 

14 

3^2 -3^4 

+ 553 

3V 4 V 77 3 

+ 156 

W 3 

+ 

256 

^2“ 77 1 

+ 231 

34 2 -4^ 3 +Tr 4 
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Equatorial True Longitudes (Cont.) 


A 

a 

A 

a 

A 

a 

THIRD SATELLITE 

-19 161 , 

V*3 

83 

l 3~ V 2 

+ 814 

'V 2 ' & 3 + 7T 4 

- 303 

2t 2 -2i 3 

+48 289 

l 3- 7T 3 

-3 736 

G 

- 333 

3^2 -3^3 

+ 47 

2 V 2 ”3 

- 114 

2 G 

19 

4 V 4 ' t 3 

+20 492 

^3^4 

- 331 

5G' -2G + 48°. 6 

8 

SVS*. 

- 217 

2-t 2 -3V"j 

-t 192 

n 3 -”4 

11 

6l 2 -6^3 

97 

U 2 -3^3 +tt 4 

+ 122 

4> -" a 

- 1 192 

V*4 

+ 150 

l 3 -2G-2n +7T 3 

-1 867 


+ 4 075 


- 1 678 

^3- 2 ^4 +7r 4 

+ 167 

^-cu 4 

+ 286 

3^-3^ 

- 100 

^2- 2Z 3 +7T l 

+ 267 

^3 “^4 

+ 67 

44 3 -« 4 

- 922 

V 2<fc 3 +7T 2 



+ 19 

s-t 3 - 5i . 

+ 1 869 

l 2~ 2l 3 +77 3 



FOURTH SATELLITE 

+ 161 

~ l 4 

- 5 775 

*4 ~ n 3 

-8 933 

G 

+ 142 


+234 697 

Z 4 ' n 4 

- 322 

2G 

- 636 

V^4 

+ 1 078 

2*4 - 2 ”4 

+ 100 

2G' -G + 173°. 17 

- 319 

2^3 - 2l 4 

+ 278 

' t 3 - 2/t 4 +7T 4 

- 781 

5G' -2G +48°. 64 

69 

3 ' t 3- 3 ' t 4 

+ 358* 

2-C -2G - 2H 

4 

+ 222 


19 

4 ' t 3- 4 ' t 4 

+ 1 764 

^ 4 -2G-2n t7T 4 

- 556 

0-^4 


* H. Andoyer [ Bull. Astr. 32, 1915 ] gives the corrected value + 386 for this coefficient. 













63 


Radii Vectors 


(projected on the equator) 

(in astronomical units of length) 


rj x 10 9 = 2.819 55 [10 6 + 

r 2 x 10 9 = 4 .486 12 [ 10 6 + 

20 cos ('Cj -^ 2 ) 

55 cos ( -i 2 -i 3 ) 

-4125 cos(2^ 1 -2^ 2 ) 

-9430 cos (it 2 - 2 ^) 

-20 cos (3^ -3^) 

-55 cos (3 ^ 2 - 3 ^ 3 ) 

-15 cos (4 , i 1 -4^ 2 ) 

-30 cos (4^ 2 -4^ 3 ) 

+ 10 cos (-ij ~'t 3 ) 

-85 cos (^ 2 -rr 2 ) 

-15 cos (2^-2^) 

-310 cos (^ 2 “ 77 3 ) 

- 45 COS (-ij -77j) 

-130 cos (^ 2 -^ 4 )] 

- 35 cos (-£ 4 -tt 3 ) 


-15 cos (-£j - 7 ^)] 


r 3 x 10 9 = 7. 155 87 [ 10 6 + 

r 4 x 10 9 = 12.586 25 [10 6 + 

635 cos(^ 2 -^ 3 ) 

100 cos(^ 3 -^ 4 ) 

+ 15 cos ( 2^ 2 — 2 ^) 

+ 20 cos ( 2^3 - 2 ^ 4 ) 

-15 cos (3 ^ 2 -3^ 3 ) 

+ 180 cos ( ^ 4 - 77 3 ) 

+ 25 cos ( i 3 --t 4 ) 

-7375 cos (^ 4 -77 4 ) 

-140 cos (2 ^ 3 -2^ 4 ) 

-25 cos (2 ^ 4 - 277 4 ) 

-10 cos ( 3^3 -3"& 4 ) 

-55 cos (^ 4 - 2 G- 2 n+ 77 4 )] 

- 1515 cos (% 3 -rr 3 ) 


- 645 cos ( ^ 3 - 7 t 4 ) ] 






Additional Terms of the Longitudes in Sampson's System 

Let Av. = 2 w x io -8 I A sin a (rad) be the correction to be added to v. in 
order to obtain the true longitude of i th satellite in Sampson's co-ordinate 
system. 


i 

A 

a 

i 

-11 667 



- 206 



78 

2^-1 p - cj 2 

2 

-11 553 

21 2 -2^j 


- 3 514 

2l 2 — 0 — co 2 


- 186 

2'C 2 -i//-a> 3 


44 

2-Cj -</>- w 4 


- 267 

2l 2 -2c^ 2 

3 

-11 025 

2^3-20 


+ 122 

2 l 3 - 4>-^ 2 


- 1 303 

2't 3 -4’-co 3 


- 264 

2*3-0-^ 

4 

-18 661 

2 ^ 4 - 2 </> 


- 1 781 

2 ' C 4 - 1 //_ oj 4 


+ 189 

2^ 4 -</'-^ 3 


83 

2 ^ 4 - 2" 4 


28 * 

2l. - 2 G - 2 tt 

4 


* Only to be introduced if H. Andoyer’s correction in the equatorial longitude is adopted. 

N. B. To obtain full-accuracy determinations of the positions, v A should be substituted to 
in the arguments a . 



65 


Latitude Series 

& . being the latitude of the i th satellite, relative to Sampson’s ( x , y) - 
plane: 

t g & i = 10“ 7 2Bsin/3 



*Erroneoosly 3-i^ — 3^ in Sampson’s Tables. 

N.B. To obtain full accuracy determinations of the positions, should be substituted to /£,. in 
the arguments / 3 ! 1 
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